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ABSTRACT 

The  impact  of  ocean  surface  waves  on  the  structure  and  intensity  of  three  tropical  cyclones  and  a  topograph¬ 
ically  forced  bora  event  is  investigated  using  a  coupled  atmosphere-ocean  wave  modeling  system.  The  coupled 
system  is  capable  of  representing  surface  momentum  fluxes  that  are  enhanced  due  to  young  ocean  waves  in 
fetch-limited  conditions,  which  yield  surface  roughness  lengths  that  significantly  depart  from  the  conventional 
Charnock-type  formulation.  In  general,  the  impact  of  ocean-wave-induced  stress  on  the  tropical  cyclone  central 
pressure  was  quite  variable  with  ocean  wave  feedback  resulting  in  changes  ranging  from  8  hPa  deeper  to  3  hPa 
shallower.  The  increased  low-level  stress  due  to  the  ocean  waves  reduced  the  near-surface  winds  by  2-3  ms-1, 
with  local  differences  in  excess  of  10  m  s”1,  which  directly  led  to  a  10%  reduction  in  the  significant  wave 
height  maxima.  The  reduced  significant  wave  heights  in  the  coupled  model  were  in  closer  agreement  with 
observations  for  Tropical  Cyclone  Bonnie  than  for  the  uncoupled  model.  The  tropical  cyclone  tracks  were 
generally  insensitive  to  ocean  wave  feedback  effects.  The  boundary  layer  structure  was  found  to  be  generally 
insensitive  to  large  roughness  enhancements  associated  with  coupled  ocean  wave  feedbacks  for  topographically 
forced  high  wind  phenomena,  such  as  the  bora. 


1.  Introduction 

The  exchange  of  heat,  momentum,  and  moisture  be¬ 
tween  the  air  and  sea  has  long  been  recognized  as  a 
fundamental  process  in  the  development  of  a  number 
of  mesoscale  atmospheric  phenomena  such  as  tropical 
cyclones,  extratropical  cyclones,  boundary  layer  jets, 
coastal  fronts,  and  precipitating  systems.  Numerical  pre¬ 
diction  of  these  phenomena  is  often  critically  dependent 
on  the  fidelity  of  surface  flux  representation  (e.g.,  Eman¬ 
uel  1995).  This  air-sea  exchange  takes  place  at  the  fluid 
interface  in  the  wave  boundary  layer,  which  underscores 
the  potential  importance  of  ocean  surface  waves  in  air- 
sea  interaction  processes. 

Ocean  surface  gravity  waves  can  have  a  significant 
impact  on  the  interaction  between  the  atmosphere  and 
ocean  through  modulation  of  the  transfer  of  momentum 
and  enthalpy  across  the  air-sea  interface  (e.g.,  Donelan 
1990).  As  air  flows  over  the  sea,  wind  stress  contributes 
to  both  the  generation  of  surface  gravity  waves  and  the 
forcing  of  local  ocean  currents  and  circulations.  The 
conventional  representation  of  the  surface  roughness  ef¬ 
fects  over  the  sea  due  to  ocean  waves,  based  on  the 
scaling  arguments  of  Charnock  (1955),  is  used  by  nearly 
all  atmospheric  research  and  operational  models,  and  is 
strictly  valid  only  for  fully  developed  ocean  wave  con- 
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ditions.  Under  high  wind  conditions,  however,  wind  di¬ 
rection  and  speed  are  often  time  dependent,  such  as  for 
a  translating  tropical  cyclone  or  topographically  driven 
flows  in  the  coastal  zone.  In  these  fetch-limited  con¬ 
ditions,  surface  ocean  waves  have  an  increasingly  im¬ 
portant  impact  on  the  momentum  flux  in  the  atmospheric 
and  oceanic  boundary  layers  (Donelan  1990).  In  situ¬ 
ations  when  the  wave  age  is  small,  the  wave-induced 
stress  composes  a  significant  fraction  of  the  total  stress. 
A  number  of  previous  studies  using  field  measurements 
have  documented  this  dependence  of  the  atmospheric 
momentum  flux  on  the  ocean  wave  age  (e.g..  Smith  et 
al.  1992).  It  should  be  noted  that  experimental  and  the¬ 
oretical  evidence  suggests  that  although  a  sea  state  de¬ 
pendence  exists  for  the  surface  momentum  flux,  no  sim¬ 
ilar  dependence  has  been  suggested  for  heat  and  mois¬ 
ture  fluxes  (Geernaert  et  al.  1987).  The  interaction  be¬ 
tween  the  sea  and  air  is  especially  complex  at  high  wind 
speeds  where  flux  exchange  processes  may  be  impacted 
by  sea  spray  and  spume  (Kepert  et  al.  1999;  Wang  et 
al.  2001)  and  airflow  separation  due  to  breaking  surface 
waves  (Kudryavtsev  and  Makin  2001). 

Wave-induced  stress  may  be  a  significant  component 
of  the  total  momentum  stress  in  the  atmospheric  bound¬ 
ary  layer  over  the  ocean  (e.g.,  Donelan  1990;  Janssen 
1991)  and  has  been  suggested  to  enhance  the  decay  of 
extratropical  systems  on  the  mesoscale  (e.g.,  Doyle 
1995;  Desjardins  et  al.  2000),  impact  the  wave  predic¬ 
tion  in  extreme  sea  states  associated  with  midlatitude 
cyclones  (Lalbeharry  et  al.  2000),  and  modulate  the  ex- 
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tratropical  storm  track  (Weber  et  al.  1993;  Janssen  and 
Viterbo  1996).  Zhang  and  Perrie  (2001)  found  that  heat 
and  moisture  fluxes  dominate  momentum  feedbacks  due 
to  wave  effects  behind  cold  fronts  in  extratropical  cy¬ 
clones,  while  momentum  fluxes  appear  to  be  more  im¬ 
portant  ahead  of  and  near  the  cyclone  center.  Powers 
and  Stoelinga  (2000)  applied  a  coupled  air-sea  model 
to  Lake  Erie  and  noted  that  roughness  enhancement  due 
to  wave  age  effects  modulated  the  local  planetary 
boundary  layer  and  downstream  flow.  The  impact  of 
wind-generated  ocean  waves  on  the  structure  and  evo¬ 
lution  of  mesoscale  phenomena  that  are  primarily  forced 
by  (i)  air-sea  fluxes  such  as  tropical  cyclones  and  (ii) 
coastal  orography  have  received  considerably  less  at¬ 
tention.  This  study  describes  a  series  of  numerical  sim¬ 
ulations  using  the  Naval  Research  Laboratory’s  (NRL) 
Coupled  Ocean-Atmosphere  Mesoscale  Prediction  Sys¬ 
tem  (COAMPS)  coupled  with  the  Wave  Model  (WAM) 
to  evaluate  the  impact  of  air-sea  interaction  on  the  de¬ 
velopment  and  track  of  tropical  cyclones  and  the  evo¬ 
lution  of  topographically  forced  flows  in  the  coastal 
zone.  Although  we  recognize  that  the  ocean  mixed  layer 
and  sea  spray  may  be  important  components  of  the  air- 
sea  interaction  processes  (e.g.,  Bao  et  al.  2000),  in  this 
study  we  focus  only  on  interaction  of  sea  surface  waves 
and  the  atmosphere  as  a  first  step  toward  fully  modeling 
the  air-sea  interaction  system.  The  numerical  model  is 
described  in  section  2.  The  results  of  the  coupled  sim¬ 
ulations  are  presented  in  section  3,  followed  by  the  sum¬ 
mary  and  conclusions  in  section  4. 


2.  Numerical  model  description 

The  atmospheric  model  used  in  this  study  is  NRL’s 
COAMPS  (Hodur  1997;  Hodur  and  Doyle  1999).  Finite- 
difference  approximations  are  used  in  the  model  to  rep¬ 
resent  the  fully  compressible,  nonhydrostatic  equations 
that  govern  atmospheric  motions.  The  equations  are 
solved  in  three  dimensions,  with  a  nested-grid  mesh 
capability,  using  a  terrain-following  vertical  coordinate 
transformation.  The  finite-difference  schemes  are  of  sec¬ 
ond-order  accuracy  in  time  and  space.  The  compressible 
equations  are  integrated  efficiently  through  the  use  of  a 
time-splitting  technique  that  features  a  semi-implicit 
treatment  in  the  vertical  for  the  acoustic  modes. 

A  prognostic  equation  for  the  turbulent  kinetic  energy 
(TKE),  e,  budget  is  used  to  parameterize  the  planetary 
boundary  layer  and  free-atmospheric  turbulent  mixing 
and  diffusion  based  on  the  level  2.5  formulation  of  Mel- 
lor  and  Yamada  (1982)  as  follows: 


de 

at 


de  de  de 

—  u - v - & - b  BP  +  SP  +  D  —  e,  (1) 

dx  dy  dcr 


where  u  and  v  are  the  horizontal  wind  velocity  com¬ 
ponents,  BP  is  the  buoyancy  production,  SP  is  the  shear 
production,  D  is  the  subgrid-scale  vertical  mixing  of 
TKE,  and  e  is  the  dissipation  rate.  The  eddy  diffusivity 
coefficients  for  heat  and  momentum  are  directly  pro¬ 


portional  to  em.  The  surface  fluxes  are  computed  fol¬ 
lowing  the  Louis  et  al.  (1982)  formulation,  which  is  an 
approximation  to  Monin-Obukhov  similarity  theory.  In 
the  standard  application  of  COAMPS  when  a  predictive 
ocean  wave  model  is  not  used  in  a  tightly  coupled  man¬ 
ner,  z0  is  represented  over  the  open  sea  by  the  Charnock 
relationship  (Charnock  1955): 

z0  =  (2) 

gP 

where  t  is  the  magnitude  of  the  turbulent  stress,  p  is 
the  density,  and  a  is  the  Charnock  parameter  taken  to 
be  0.0185. 

A  force-restore  method  is  used  to  parameterize  the 
land  surface  energy  budget.  Parameterization  of  short- 
and  longwave  radiation  processes,  and  subgrid-scale 
moist  convective  processes,  are  included  as  well.  Budget 
equations  for  cloud  water,  cloud  ice,  raindrops,  snow¬ 
flakes,  and  water  vapor  are  used  to  represent  the  grid- 
scale  evolution  of  the  moist  processes.  For  more  details 
of  the  parameterization  methodology,  see  Hodur  (1997) 
and  Hodur  and  Doyle  (1999). 

Multivariate  optimum  interpolation  analysis  of  upper- 
air  sounding,  surface,  aircraft  (automated  weather  re¬ 
ports  from  commercial  aircraft;  ACARS)  and  satellite 
data  that  are  quality  controlled  is  the  basis  for  the  initial 
state  of  the  nonhydrostatic  model.  The  Navy  Opera¬ 
tional  Global  Analysis  and  Prediction  System  (NO¬ 
GAPS)  forecast  fields  are  used  for  lateral  boundary  con¬ 
ditions.  The  application  of  COAMPS  in  this  type  of 
manner  closely  emulates  a  real-time  numerical  predic¬ 
tion  system  in  spite  of  the  simulations  being  performed 
in  a  hindcast  mode.  The  atmospheric  topographic  data 
are  based  on  the  U.S.  Defense  Mapping  Agency’s  100- 
m  resolution  dataset. 

The  atmospheric  model  is  coupled  in  an  interactive 
mode  to  the  third  generation  (cycle  4)  of  the  WAM 
(WAMDI  Group  1988).  The  quasi-linear  theory  of 
wind- wave  generation  following  Janssen  (1989)  and 
Janssen  (1991)  is  applied.  The  transport  equation  is  in¬ 
tegrated  to  solve  for  the  wave  variance  spectrum  for  25 
frequencies  and  24  directions.  The  processes  repre¬ 
sented  in  WAM  include  wind-wave  interaction,  transfer 
of  energy  within  the  wave  spectrum,  wave  propagation, 
and  dissipation  by  whitecapping. 

The  coupling  methodology  follows  Janssen  (1989) 
and  Janssen  (1991)  and  includes  the  processes  repre¬ 
sented  by  mutual  interaction  of  the  wind  waves  and 
boundary  layer  stress.  The  roughness  length  used  in  the 
coupled  simulation  is  represented  by 


T 

P gp(  1  -  tJt)0S 


(3) 


where  t  is  the  total  stress  and  tw  is  the  wave-induced 
stress.  The  constant  j 8  is  chosen  as  0.01  implying  that 
(3)  reduces  to  the  standard  Charnock  relationship  (2) 
for  a  saturated  wave  state  (rw  =  0).  The  wave-induced 
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stress  is  primarily  determined  by  the  high  frequency 
portion  of  the  spectrum  because  these  waves  have  the 
fastest  growth  rate  due  to  the  wind  forcing  (Janssen 
1989).  For  a  young  wind  sea,  it  follows  that  the  effective 
Charnock  parameter  can  be  enhanced  by  an  order  of 
magnitude.  The  simultaneous  coupling  is  physically 
achieved  through  communicating  the  10-m  wind  speed 
to  the  wave  model  every  WAM  time  step.  An  iterative 
technique  is  then  used  to  calculate  t„  based  on  the  wind 
speed,  drag  coefficient,  and  stress.  The  new  wave-in¬ 
duced  Charnock  parameter  is  subsequently  used  in  the 
Louis  et  al.  (1982)  flux  calculations  in  the  atmospheric 
model.  This  coupling  technique  using  the  WAM  and  an 
atmospheric  model  has  been  applied  previously  to  study 
air-sea  interaction  in  extratropical  cyclones  (e.g.,  Doyle 
1995;  Lionello  et  al.  1998;  Desjardins  et  al.  2000;  Lal- 
beharry  et  al.  2000)  and  tropical  cyclones  (Bao  et  al. 
2000).  The  coupling  technique  using  the  WAM  is  cur¬ 
rently  being  applied  operationally  at  the  European  Cen¬ 
tre  for  Medium-Range  Weather  Forecasts  (ECMWF) 
and  the  results  indicate  a  substantial  positive  impact  on 
the  skill  of  the  atmosphere  and  ocean  wave  model  fore¬ 
casts  (Janssen  et  al.  2002).  Desjardins  et  al.  (2000) 
showed  that  an  empirical  sea-state-dependent  z0,  based 
on  data  collected  during  the  Humidity  Exchange  Over 
the  Sea  (HEXOS)  experiment,  yields  consistent  results 
with  (3).  In  this  study,  the  uncoupled  simulations  are 
performed  with  simultaneous  integration  of  the  WAM 
except  the  standard  Charnock  parameter  of  0.0185  is 
used  in  the  atmospheric  model  flux  calculations. 

3.  Coupled  simulations 

To  isolate  the  importance  of  the  wave  coupling,  two 
experiments  are  performed  for  each  case.  The  first  sim¬ 
ulation  makes  use  of  the  conventional  Charnock  param¬ 
eter  in  the  atmosphere  and  includes  a  wave  simulation 
based  on  the  atmospheric  winds  every  WAM  time  step. 
The  second  simulation  includes  the  wave-induced  stress 
effects  in  the  modified  Charnock  formulation  following 
(3)  and  also  includes  the  feedback  of  the  atmospheric 
winds,  already  modified  by  wave-induced  roughness, 
back  to  the  waves. 

Coupled  and  uncoupled  COAMPS/WAM  simulations 
were  performed  for  Tropical  Cyclones  Bonnie,  Mitch, 
and  Bret,  and  a  topographically  forced  high  wind  event. 
A  horizontal  grid  dimension  of  133  X  133  with  an  in¬ 
crement  of  12  km  was  used  for  both  CO  AMPS  and 
WAM  model  simulations  of  Tropical  Cyclones  Mitch 
and  Bret.  Two  nested  grids  were  used  for  the  simulations 
of  Bonnie.  A  fine  mesh  grid  of  181  X  181  points  and 
an  increment  of  6  km  were  used  with  an  identical  grid 
domain  and  resolution  for  the  single  WAM  mesh.  The 
coarse  mesh  for  the  atmospheric  model  contained  133 
X  133  points  with  a  resolution  of  18  km.  The  simulation 
of  the  topographically  forced  wind  event  made  use  of 
a  single  mesh  with  133  X  133  points  and  a  horizontal 
grid  increment  of  5  km.  Thirty  vertical  levels  were  used 


in  the  atmospheric  model  with  a  top  at  34  km  and  a 
lowest  model  level  at  10  m  above  ground.  A  terrain¬ 
following  vertical  coordinate  is  used  that  features  hor¬ 
izontal  model  surfaces  over  water  and  flat  land  (Hodur 
1997).  In  the  WAM,  the  wave  spectrum  is  discretized 
using  24  directions  and  25  frequencies  in  the  range 
0.033-0.33  Hz  for  all  experiments.  The  directional  dis¬ 
cretization  is  rotated  by  1.5°  to  avoid  a  systematic  bias 
in  the  energy  propagation. 

a.  Tropical  Cyclone  Bonnie  simulations 

Tropical  Cyclone  Bonnie  strengthened  to  hurricane 
force  at  0600  UTC  22  August  1998  when  it  was  located 
400  km  north  of  Hispaniola.  Bonnie  moved  on  a  general 
west-northwest  heading  reaching  maximum  strength  of 
51  m  s  1  and  a  central  pressure  of  954  hPa.  Eventually, 
Bonnie  made  landfall  near  Wilmington,  North  Carolina, 
at  approximately  0330  UTC  27  August.  The  time  period 
of  interest  in  this  study  is  during  an  earlier  phase  of  the 
tropical  cyclone  on  24  August.  On  this  day,  the  sea 
surface  directional  wave  spectrum  was  measured  using 
the  National  Aeronautics  and  Space  Administration 
(NASA)  airborne  Scanning  Radar  Altimeter  (SRA)  car¬ 
ried  aboard  a  National  Oceanic  and  Atmospheric  Ad¬ 
ministration  (NOAA)  WP-3D  hurricane  research  aircraft 
at  a  height  of  1.5  km  (Wright  et  al.  2001).  A  description 
of  the  instrument,  error  characteristics,  and  applications 
can  be  found  in  Walsh  et  al.  (1996)  and  Wright  et  al. 
(2001). 

The  Charnock  parameter  for  1200  UTC  24  August 
1998  (24-h  simulation  time)  for  the  coupled  simulation, 
shown  in  Fig.  la,  indicates  a  large  region  of  enhanced 
Charnock  parameter  outside  the  tropical  cyclone  eye. 
The  northern  or  right-of-track  half  of  the  tropical  cy¬ 
clone  contains  the  region  with  the  largest  wave-induced 
stress,  resulting  in  Charnock  parameter  maxima  in  ex¬ 
cess  of  0.09.  The  wave  age  parameter,  Cp/u *  where  Cp 
is  the  phase  speed  of  the  peak  spectral  frequency  and 
m*  is  the  friction  velocity,  for  1200  UTC  24  August,  is 
shown  in  Fig.  lb.  The  wave  age  is  consistent  with  the 
Charnock  parameter  distribution  and  indicates  that  the 
youngest  ocean  waves  are  located  in  the  northern  half, 
with  swell  dominating  the  eye  and  outer  regions  of  the 
tropical  cyclone.  Clearly  in  these  tropical  cyclone  sit¬ 
uations,  the  strong  curvature  of  the  wind  field  limits  the 
fetch.  As  a  result  of  the  fetch-limited  conditions,  the 
roughness  length  was  enhanced  substantially  from  val¬ 
ues  of  10-5  m  present  in  the  near-calm  regions  within 
the  eye  and  outside  of  the  tropical  cyclone,  to  greater 
than  0.03  m  in  the  northeastern  quadrant  (Fig.  lc). 

Typically,  the  wind  speed  in  a  tropical  cyclone  is 
largest  in  the  region  that  contains  winds  with  the  same 
direction  as  the  storm  motion  because  the  forward  trans¬ 
lation  speed  of  the  storm  contributes  to  higher  speeds. 
In  this  region,  the  waves  propagate  in  the  direction  of 
the  storm  motion  so  that  a  given  wave  packet  may  be 
under  the  influence  of  a  consistent  wind  direction  and 
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Fig.  1.  Coupled  simulation  for  1200  UTC  24  Aug  1998  of  (a) 
nondimensional  Charnock  parameter  (X10~2),  (b)  wave  age,  and  (c) 
roughness  length  (X10-2  m).  Winds  at  10  m  are  indicated  by  the 
vectors. 


Fig.  2.  Simulated  dimensionless  roughness  length  (z0*  =  gzjul ) 
as  a  function  of  the  wave  age  for  the  coupled  simulation  at  1200 
UTC  24  Aug.  The  Charnock  parameter  of  0.0185  is  shown  by  the 
dashed  line  at  ln(z0*)  —  —  4.0.  The  solid  line  represents  the  Smith  et 
al.  (1992)  best  fit. 


speed  for  a  longer  period  and  distance  (Wright  et  al. 
2001).  These  resonance  processes  act  to  increase  the 
fetch  and  duration  in  the  same  direction  as  the  storm 
motion  (Young  1988;  MacAfee  and  Bowyer  2000),  par¬ 
ticularly  when  the  group  velocity  of  the  waves  and  trans¬ 
lation  speed  of  the  tropical  cyclone  are  comparable 
(Moon  et  al.  2002).  These  factors  along  with  the  fetch- 
limited  effects  due  to  curvature  in  the  wind  field  result 
in  an  enhanced  wave-induced  surface  roughness  in  the 
right  half  of  the  tropical  cyclone.  Additionally,  the 
asymmetry  in  the  storm-scale  structure,  with  the  greatest 
precipitation  and  cloud  shield  located  in  the  northern 
half,  as  indicated  by  satellite  imagery,  may  also  have 
influenced  the  asymmetric  wave  state  characteristics. 

The  relationship  between  the  simulated  wave  age  pa¬ 
rameter  and  the  dimensionless  roughness  length  ( gz0 / 
u%)  for  1200  UTC  24  August  is  shown  in  Fig.  2.  Model- 
simulated  data  extracted  from  every  fifth  grid  point  over 
the  entire  domain  are  displayed.  The  regression  fit  based 
on  observational  data  from  Smith  et  al.  (1992)  taken 
during  HEXOS,  shown  in  Fig.  2,  indicates  a  reasonable 
agreement  between  the  observations  and  simulated 
wave  age  and  Charnock  parameters,  as  is  also  discussed 
in  the  coupled  model  results  of  Janssen  et  al.  (2002). 
In  this  case  of  extreme  winds  the  young  waves  in  the 
range  from  approximately  2.3  to  3.0  appear  to  have  a 
larger  negative  slope  than  the  Smith  et  al.  (1992)  data 
would  indicate.  The  larger  negative  slope  exhibited  by 
this  distribution  is  in  agreement  with  the  wave-age-de- 
pendent  Charnock  parameters  computed  by  Oost  et  al. 
(2001)  based  on  eddy  correlation  measurements  in  the 
North  Sea.  The  uncoupled  simulations,  which  do  not 
feedback  the  wave-induced  stress,  use  a  dimensionless 
roughness  length  of  —4.0  or  a  =  0.0185  over  the  water. 
A  constant  Charnock  parameter  of  0.0185  is  sufficiently 
accurate  for  a  fully  saturated  wave  state  ( Cplu *  >  25). 
However,  for  a  relatively  young  wind  sea  ( Cplu *  <  25), 
the  use  of  the  standard  Charnock  value  significantly 
underestimates  the  dimensionless  roughness  (e.g..  Smith 
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Fig.  3.  Neutral  drag  coefficient  (X10-3)  as  a  function  of  wind 
speed  from  the  coupled  simulation  at  1200  UTC  24  Aug.  Gray-shaded 
lines  indicate  drag  coefficients  associated  with  very  young  and  fully 
developed  wave  states  following  Smith  et  al.  (1992). 

et  al.  1992).  It  follows  that  young  waves  will  be  as¬ 
sociated  with  a  larger  effective  drag  coefficient  than  for 
a  fully  developed  wave  state  (e.g.,  Geernaert  et  al. 
1987),  as  illustrated  in  Fig.  3.  Thus,  the  differences 
between  the  uncoupled  and  coupled  simulations  can  be 
attributed  to  the  roughness  enhancement  due  to  young 
ocean  waves.  It  should  be  noted  that  in  this  simulation 
the  dimensionless  roughness  exceeds  the  Charnock  val¬ 
ue  for  fully  developed  waves  at  a  number  of  points  in 
general  agreement  with  the  Oost  et  al.  (2001)  North  Sea 
measurements.  This  may  not  be  unexpected  given  the 
complex  wave  state  that  dominates  nearly  the  entire 
inner-mesh  domain.  Applications  of  the  coupled  model 
using  this  parameterization  to  other  cases  such  as  ex- 
tratropical  cyclones,  which  contain  weaker  wind  speeds, 
indicate  a  closer  agreement  with  the  Smith  et  al.  (1992) 
regression  curve  for  the  fully  developed  regime  (e.g., 
see  Fig.  4  in  Doyle  1995). 

Relevant  aspects  of  the  Bonnie  simulation  for  1200 
UTC  24  August  (24-h  simulation  time)  and  0000  25 
August  (36-h  simulation  time)  are  summarized  in  Table 
1.  The  central  pressures  for  the  coupled  and  uncoupled 
simulations  are  quite  similar  and  within  2  hPa  of  each 
other,  a  result  also  noted  by  Desjardins  et  al.  (2000)  for 
extratropical  cyclones.  However,  the  observed  central 
pressure  is  15-20  hPa  deeper  than  simulated  in  the  cou¬ 
pled  and  uncoupled  simulations.  The  tropical  cyclone 


track  is  nearly  identical  between  the  coupled  and  un¬ 
coupled  simulations.  Similar  maximum  horizontal  ve¬ 
locities  at  10  m  and  500  m  are  attained  in  the  two  sim¬ 
ulations. 

The  most  substantial  differences  between  the  simu¬ 
lations  are  apparent  in  the  significant  wave  height  fields. 
The  coupled  simulation  has  a  significant  wave  height 
maximum  that  is  reduced  by  approximately  15%  relative 
to  the  uncoupled  case.  Both  the  coupled  and  uncoupled 
simulations  make  use  of  new  atmospheric  winds  every 
wave  model  time  step.  The  SRA-derived  significant 
wave  height  and  dominant  wavelength  superimposed  on 
the  NOAA’s  Atlantic  Oceanographic  and  Meteorological 
Laboratory/Hurricane  Research  Division  (AOML/ 
HRD)  surface  wind  field  analysis  are  shown  in  Fig.  4 
(from  Wright  2000).  The  maximum  significant  wave 
height  derived  from  this  composite  analysis  of  five  flight 
segments  is  approximately  10.5  m.  It  is  noteworthy  that 
the  dominant  waves  propagate  at  significant  angles  to 
the  low-level  wind  direction.  The  significant  wave 
height  maximum  is  located  in  the  northeast  quadrant 
(Fig.  4a)  where  the  field  is  dominated  by  swell  (Fig. 
4b).  In  the  southern  half  of  the  tropical  cyclone,  the 
dominant  wavelength  is  considerably  shorter  relative  to 
the  long  wavelength  associated  with  the  swell  located 
in  the  northern  half. 

The  simulated  significant  wave  heights  for  the  un¬ 
coupled  and  coupled  simulations  of  Tropical  Cyclone 
Bonnie  are  shown  in  Fig.  5.  As  noted  above,  the  sig¬ 
nificant  wave  height  maximum  in  the  northeastern  quad¬ 
rant  is  considerably  reduced  in  the  coupled  simulation, 
in  spite  of  the  dominance  of  the  swell.  The  coupled 
simulation  is  in  overall  closer  agreement  with  the  SRA- 
derived  wave  field  than  the  uncoupled  case.  Lalbeharry 
et  al.  (2000)  also  found  that  two-way  coupling  has  some 
positive  impact  through  the  mitigation  of  the  overpred¬ 
iction  of  significant  wave  height  maxima  apparent  in 
uncoupled  simulations,  particularly  in  cases  with  ex¬ 
treme  sea  states  as  considered  here.  The  dominant  wave¬ 
length  is  shown  in  Fig.  6  for  the  uncoupled  and  coupled 
simulations.  Considering  the  complexity  of  both  the 
wave  and  wind  fields,  the  simulated  wavelength  contains 
many  common  characteristics  with  the  SRA  composite 
wavelength  field  (Fig.  4b)  including  the  dominance  of 


Table  1.  Summary  of  the  simulated  central  pressure  (SLP),  wind  speed  maximum  (Umax),  significant  wave  height  maximum  (Hs),  and 
maximum  Charnock  parameter  value  (amax)  for  Tropical  Cyclone  Bonnie  valid  at  1200  UTC  24  Aug  (24  h)  and  0000  UTC  25  Aug  1998 
(36  h).  Observed  estimates  from  the  best  track  are  shown  as  well.  The  uncoupled  and  coupled  wind  speed  maxima  are  shown  for  the  10- 
and  500-m  levels,  respectively. 


Variable 

Simulation  time  (h) 

Uncoupled 

Coupled 

Observed 

SLP 

24 

982.6  hPa 

983.2  hPa 

962  hPa 

SLP 

36 

980.2  hPa 

978.7  hPa 

964  hPa 

24 

35.4/54.1  m  s"1 

32.1/52.2  m  s”1 

51.0  m  s-1 

u, „ 

36 

33.4/49.7  m  s“‘ 

35.1/54.7  m  s'1 

51.0  m  s~‘ 

H, 

24 

12.8  m 

10.9  m 

10.5  m 

Hs 

36 

12.6  m 

11.0  m 

24 

0.0185 

0.15 

amax 

36 

0.0185 

0.11 
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Fig.  4.  NASA  airborne  SRA-derived  (a)  significant  wave  height  (m)  and  (b)  dominant  wavelength  (m)  based  on 
NOAA  WP-3D  research  aircraft  data  from  24  Aug  1998  (modified  from  Wright  2000  and  Wright  et  al.  2001). 


swell  in  the  northern  half  and  shorter  wavelengths  in 
the  southern  portion  of  the  tropical  cyclone.  The  dom¬ 
inant  wavelength  in  the  coupled  simulation  (Fig.  6b)  is 
slightly  shorter  in  the  northern  half  of  the  tropical  cy¬ 
clone  than  the  uncoupled  simulation  (Fig.  6a).  Both  the 
coupled  and  uncoupled  simulations  contain  a  region 
south  of  the  eye  characterized  by  wavelengths  less  than 
100  m,  whereas  the  SRA  observations  indicate  slightly 
longer  wavelengths  greater  than  150  m.  The  SRA  mea¬ 
surements  indicate  that  the  swell  propagated  in  a  north¬ 
westerly  direction  aligned  with  the  track  of  the  tropical 
cyclone.  Wright  et  al.  (2001)  hypothesized  that  the  dom¬ 
inance  of  the  northwesterly  swell  was  due  to  higher  wind 
speeds  and  wave  resonance  effects  that  increased  the 
effective  fetch,  which  occurred  in  the  northern  half  of 
the  tropical  cyclone.  The  coupled  and  uncoupled  sim¬ 
ulations  both  capture  the  many  characteristics  of  the 
SRA  directional  wave  spectra,  such  as  the  dominance 
of  northwesterly  propagating  swell  and  the  tendency  for 


multiple  modes  in  the  southeastern  quadrant.  It  is  note¬ 
worthy,  however,  that  the  radial  distance  from  the  trop¬ 
ical  cyclone  eye  to  the  wave  height  maximum  was  ob¬ 
served  to  be  approximately  70  km  based  on  the  SRA 
data  and  200  km  in  the  model  simulation.  Although  the 
coupled  simulation  appears  to  be  in  general  agreement 
with  a  number  of  observed  features  and  characteristics, 
fundamental  deficiencies  in  the  model  simulation  of  the 
tropical  storm  structure  and  intensity  clearly  exist. 

The  directional  wave  spectra  as  measured  from  the 
SRA,  shown  in  Fig.  7  (from  Wright  et  al.  2001),  illus¬ 
trate  a  number  of  complex  characteristics  of  the  wave 
field.  At  several  locations  wave  spectra  were  extracted 
from  the  coupled  simulation  for  similar  regions  as  the 
SRA  observations  (Fig.  7a)  and  are  shown  in  Fig.  8  (in 
eye-relative  coordinates).  In  the  northwestern  quadrant 
at  location  point  1,  the  observed  (Fig.  7b)  and  simulated 
(Fig.  8a)  directional  wave  spectra  contain  common  char¬ 
acteristics  including  a  westward  propagating  swell  with 


Fig.  5.  Simulated  significant  wave  height  (m)  and  direction  (vectors)  for  the  (a)  uncoupled 
and  (b)  coupled  simulations  valid  at  1200  UTC  24  Aug  1998  (24  h). 
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Fig.  6.  Wavelength  (m)  of  the  wave  spectrum  peak  valid  at  1200  UTC  24  Aug  1998  (24  h)  for 
the  (a)  uncoupled  and  (b)  coupled  models. 
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Fig.  7.  Observed  directional  wavenumber  spectra  from  the  NASA  SRA  measured  between  2030  and  2115  UTC  24 
Aug  1998.  (top-left)  The  NO  A  A  aircraft  track  and  HRD  surface  winds  field  (speed  of  30  m  s"1  corresponds  to  a 
length  of  10  km)  and  the  location  of  the  spectra  of  interest  relative  to  the  tropical  cyclone  center  (modified  from 
Wright  et  al.  2001).  The  outer  solid  circle  in  each  spectrum  corresponds  to  100  m,  and  the  inner  circles  represent 
wavelengths  of  150,  200,  250,  300,  and  350  m. 
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Fig.  8.  Simulated  directional  wave  spectra  from  the  coupled  simulation  valid  at  1200  UTC  24  August  1998.  The 
solid  outer  circle  corresponds  to  the  50-m  wavelength.  The  dashed  circles  indicate  the  100-,  150-,  200-,  250-,  300-, 
and  350-m  wavelengths.  The  spectra  are  individually  normalized.  The  locations  correspond  to  the  relative  positions 
to  the  tropical  cyclone  as  shown  in  Fig.  7a  and  are  denoted  by  the  number  in  the  lower  right. 


wind-generated  waves  manifested  as  a  tail  on  the  spec¬ 
trum  that  extends  toward  shorter  wavelengths  or  higher 
frequencies.  The  model-simulated  peak  in  the  energy 
spectrum  occurs  near  275  m,  which  is  a  slightly  longer 
wavelength  than  the  observed  peak  at  250  m.  Near  the 
tropical  cyclone  eye  at  location  point  2,  the  observed 
(Fig.  7c)  and  simulated  (Fig.  8b)  spectra  both  contain 
a  similar  signature  with  westward  propagating  swell. 
The  wavelength  at  the  peak  of  the  spectrum,  275  m,  is 
once  again  slightly  longer  than  the  derived  value  from 
the  SRA,  250  m.  In  the  southeastern  quadrant  at  location 
point  3,  the  observed  (Fig.  7d)  and  simulated  (Fig.  8c) 
spectra  indicate  a  bimodal  wave  field  with  a  northeast 
propagating  wind-generated  sea  superimposed  on  a 
northwest  propagating  swell.  However,  at  location  3,  the 
SRA-observed  wave  state  is  characterized  by  greater 
energy  and  longer  wavelengths  in  the  northeast  prop¬ 


agating  waves  than  the  coupled  simulation  indicated. 
This  difference  may  arise  from  the  deficiencies  in  the 
atmospheric  model  simulation  of  the  tropical  cyclone 
structure.  An  additional  point  extracted  south  of  the  cy¬ 
clone  eye,  at  location  4,  indicates  an  even  more  complex 
wave  state  with  multiple  wind  sea  modes  coincident 
with  a  northwest  propagating  swell.  Although  the  iden¬ 
tical  frequency  spectrum  was  not  sampled  by  the  SRA 
at  this  point,  a  number  of  complex  trimodal  wave  spectra 
were  noted  in  the  vicinity  by  Wright  et  al.  (2001).  Over¬ 
all,  wave  spectra  from  the  uncoupled  simulation  tended 
to  contain  qualitatively  similar  structures  as  the  fully 
coupled  model  with  much  greater  differences  in  the 
wave  spectra  between  the  quadrants  surrounding  the 
tropical  cyclone  than  the  differences  between  the  two 
simulations. 

Although  significant  quantitative  differences  exist  in 
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Fig.  9.  Radial  distribution  of  the  azimuthally  averaged  sensible  and 
latent  heat  fluxes  for  the  simulation  of  Tropical  Cyclone  Bonnie  valid 
at  1200  UTC  24  Aug  (24  h).  The  uncoupled  and  coupled  simulation 
results  are  denoted  by  the  dark  and  gray  curves,  respectively. 

the  wave  properties  between  the  uncoupled  and  coupled 
simulations,  the  surface  heat  fluxes  are  quite  similar 
between  the  two  experiments.  The  radial  distribution  of 
the  sensible  and  latent  heat  flux  averaged  azimuthally 
at  the  surface  is  shown  in  Fig.  9  for  1200  UTC  24 
August  1998  (24  h).  The  averaged  latent  heat  flux  in 
the  coupled  simulation  is  approximately  10%-20%  larg¬ 
er  than  the  uncoupled  experiment  at  radial  distances  of 
50-100  km  from  the  cyclone  center.  The  enhanced  sur¬ 
face  moisture  heat  flux  is  a  direct  result  of  the  increased 
stress  associated  with  young  ocean  waves  and  the  cor¬ 
responding  augmented  Charnock  parameter.  However, 
surface  flux  parameterizations  commonly  applied  in  at¬ 
mospheric  models  such  as  COAMPS  are  not  designed 
for  the  high  wind  speed  regime  and  it  follows  that  con¬ 
siderable  uncertainties  exist  in  the  air-sea  flux  repre¬ 
sentation. 

b.  Tropical  Cyclone  Mitch  simulations 

The  observed  minimum  central  pressure  in  Tropical 
Cyclone  Mitch  was  estimated  to  be  905  hPa  with  sus¬ 
tained  surface  wind  speed  in  excess  of  80  m  s_1,  dis¬ 
tinguishing  Mitch  as  the  strongest  October  hurricane  in 
the  Atlantic  basin  on  record.  Tropical  Cyclone  Mitch 
was  responsible  for  over  9000  deaths,  making  it  one  of 
the  deadliest  Atlantic  tropical  cyclones  in  history. 


The  simulated  central  pressure,  10-  and  500-m  wind 
speed  maxima,  maximum  significant  wave  height,  and 
maximum  Charnock  parameter  for  0000  UTC  26  Oc¬ 
tober  (24-h  simulation  time)  and  0000  UTC  27  October 
(48 -h  simulation  time)  1998  are  summarized  in  Table 
2.  The  central  pressure  is  approximately  8  hPa  lower  in 
the  coupled  simulation  than  in  the  uncoupled  case  at  the 
48-h  time.  As  a  result  of  the  enhanced  surface  stress, 
the  momentum  flux  increases  considerably  due  to  the 
coupled  response.  In  the  parameterization  of  the  surface 
fluxes,  the  heat  and  moisture  fluxes  are  influenced  by 
the  wave  state  through  the  surface  momentum  flux.  For 
example,  at  0000  UTC  27  October  wave-induced  effects 
increase  the  maximum  surface  sensible  heat  flux  by 
nearly  a  factor  of  2,  from  a  maximum  of  506  to  946  W 
m  2 .  However,  uncertainties  clearly  exist  in  the  surface 
flux  parameterization  in  this  environment.  The  azi¬ 
muthally  averaged  latent  heat  flux  nearly  doubles  from 
1500  W  m~2  in  the  uncoupled  experiment  to  2900  W 
m~2  in  the  coupled  simulation  at  a  radial  distance  of  24 
km.  This  increase  in  the  heat  flux  contributes  to  a  30% 
increase  in  the  azimuthally  averaged  rainfall  that  leads 
to  an  8-hPa  reduction  in  central  pressure  in  the  coupled 
simulation.  However,  because  of  the  increased  drag  due 
to  the  young  ocean  waves,  the  maximum  wind  speed 
below  500  m  is  less  for  the  coupled  simulation.  A  second 
coupled  simulation  was  performed  using  the  surface 
heat  and  moisture  flux  computed  using  a  surface  stress 
based  on  the  Charnock  parameter  and  independent  of 
wave  age.  The  simulated  minimum  pressure  was  941.4 
hPa  for  this  coupled  simulation,  which  is  8  hPa  higher 
than  the  uncoupled  simulation  and  16  hPa  greater  than 
the  coupled  simulation  with  the  heat  and  moisture  fluxes 
dependent  on  the  wave  state.  Clearly,  the  heat  and  mois¬ 
ture  exchange  coefficients  appear  to  have  a  very  im¬ 
portant  effect  on  tropical  cyclone  intensity  relative  to 
momentum  flux  changes  due  to  wave-induced  stress,  a 
similar  conclusion  to  that  noted  by  Zhang  and  Perrie 
(2001)  for  an  extratropical  cyclone. 

As  a  result  of  weaker  surface  winds  due  to  enhanced 
roughness  associated  with  the  young  ocean  waves,  the 
significant  wave  height  maximum  is  2  m  less  in  the 
coupled  simulation,  as  shown  in  Table  2  and  Fig.  10  for 
0000  UTC  27  October.  The  increased  stress  associated 
with  the  young  ocean  waves  does  not  appear  to  have  a 


Table  2.  Same  as  Table  1,  except  for  Tropical  Cyclone  Mitch  valid  at  0000  UTC  26  Oct  (24  h)  and  0000  UTC  27  Oct  1998  (48  h). 
Observed  estimates  from  the  best  track  are  shown  as  well. 


Variable 

Simulation  time  (h) 

Uncoupled 

Coupled 

Observed 

SLP 

24 

960.1  hPa 

958.2  hPa 

923  hPa 

SLP 

48 

933.6  hPa 

925.7  hPa 

910  hPa 

U, m 

24 

41.1/62.6  m  s~ 1 

36.6/59.9  m  s'1 

67.0  m  s-1 

u, „ 

48 

52.6/77.0  m  s“‘ 

50.1/73.8  m  s”1 

77.0  m  s-‘ 

H, 

24 

14.9  m 

13.9  m 

Hs 

48 

17.2  m 

15.1  m 

24 

0.0185 
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Fig.  10.  Significant  wave  height  (m)  and  direction  (vectors)  for  (a)  uncoupled  and  (b)  coupled  simulations  of 
Tropical  Cyclone  Mitch  valid  at  0000  UTC  27  Oct  1998  (48-h  time). 


large  impact  on  the  spatial  character  of  the  ocean  wave 
response.  Unfortunately,  no  observed  estimates  of  the 
wave  fields  are  available  for  this  time.  The  Charnock 
parameter  is  shown  in  Fig.  1 1  for  0000  UTC  26  October. 
The  maximum  Charnock  parameter  is  0.1  at  this  time, 
which  is  nearly  an  order  of  magnitude  larger  than  the 
conventional  parameter  valid  for  the  open  ocean  fea¬ 
turing  a  mature  wave  state.  Also,  the  Charnock  param¬ 
eter  maximum  undergoes  a  threefold  increase  between 
0000  UTC  26  October  and  0000  UTC  27  October  co¬ 
inciding  with  the  time  period  of  rapid  intensification  of 
the  tropical  cyclone  (i.e.,  minimum  sea  level  pressure 
decreases  by  27  hPa  in  the  33-h  period  in  the  coupled 
simulation).  It  is  noteworthy  that  the  Charnock  param¬ 
eter  is  a  maximum  near  the  eyewall  structure  coincident 


Fig.  11.  Charnock  parameter  (X102)  for  Tropical  Cyclone  Mitch 
valid  at  0000  UTC  26  Oct  1998  (24  h)  shown  by  the  gray  shading 
scale.  The  0.0181  contour  is  also  shown. 


with  the  highest  wind  speeds  and  a  minimum  within  the 
eye.  Similar  to  the  simulation  of  Tropical  Cyclone  Bon¬ 
nie,  the  enhanced  roughness  effects  due  to  the  young 
wave  state  have  a  minimal  impact  on  the  tropical  cy¬ 
clone  track,  in  spite  of  the  significant  changes  in  the 
central  pressure.  For  example,  at  0000  UTC  27  October 
(48  h),  the  position  of  the  tropical  cyclone  differs  by 
one  grid  point  (12  km)  in  the  two  simulations. 

c.  Tropical  Cyclone  Bret  simulations 

Tropical  Cyclone  Bret  was  an  intense  hurricane  with 
a  relatively  small  size  that  formed  in  the  Gulf  of  Mexico 
and  eventually  made  landfall  along  the  south  Texas 
coast.  Near  the  time  of  peak  intensity,  hurricane  force 
winds  were  confined  to  a  narrow  radius  of  approxi¬ 
mately  50  km  wide  in  the  northern  portion  of  the  cyclone 
and  25  km  to  the  south. 

The  simulated  characteristics  of  the  atmospheric  and 
wave  fields  are  summarized  in  Table  3.  The  central  pres¬ 
sure  of  the  tropical  cyclone  is  approximately  3  hPa  lower 
in  the  coupled  simulation  than  in  the  uncoupled  case, 
in  slightly  better  agreement  with  the  estimated  values 
based  on  in  situ  observations.  The  maximum  wind  speed 
at  10  and  500  m  is  once  again  slightly  reduced  due  to 
the  wave  age  effects.  The  significant  wave  height  max¬ 
imum  is  approximately  0.5  m  lower  in  the  coupled  sim¬ 
ulation,  generally  similar  to  the  response  of  Tropical 
Cyclones  Bonnie  and  Mitch.  One  noteworthy  aspect  of 
Bret  was  that  the  operational  NOGAPS  had  a  large  track 
error  of  over  400  km  at  48  h,  verifying  at  1200  UTC 
22  August  1999.  The  uncoupled  and  coupled  simula¬ 
tions  contained  a  similar  error,  implying  that  the  track 
simulation  is  once  again  insensitive  to  any  wave  cou¬ 
pling  issues. 
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Table  3.  Same  as  Table  1,  except  for  Tropical  Cyclone  Bret  valid  at  1200  UTC  21  Aug  (24  h)  and  1200  UTC  22  Aug  1999  (48  h). 


Variable 

Simulation  time  (h) 

Uncoupled 

Coupled 

Observed 

SLP 

24 

987.2  hPa 

983.9  hPa 

979  hPa 

SLP 

48 

955.2  hPa 

952.1  hPa 

944  hPa 

£/, „ 

24 

30.4/45.7  m  s"1 

30.5/47.4  m  s'1 

41.0  m  s-1 

u ^ 

48 

44.1/66.9  m  s-1 

41.8/66.2  m  s“> 

64.0  m  s-1 

H, 

24 

8.3  m 

7.6  m 

H, 

48 

12.3  m 

11.8  m 

24 

0.0185 

0.10 

amax 

48 

0.0185 

0.32 

d.  Bora  simulations 

The  Mesoscale  Alpine  Programme  (MAP)  took  place 
during  fall  1999  with  the  objective  to  measure  and  un¬ 
derstand  topographically  forced  moist  and  dry  atmo¬ 
spheric  phenomena  in  the  vicinity  of  the  Alps  (Bou- 
geault  et  al.  2001).  On  7  November  1999,  a  severe 
downslope  wind  storm  occurred  in  the  lee  of  the  Dinirac 
Alps  of  Slovenia  and  Croatia;  this  type  of  wind  often 
referred  to  locally  as  a  bora  (Smith  1987).  Research 
aircraft  measurements  in  the  low  levels  indicate  wind 
speeds  were  in  excess  of  30  m  s_1  along  the  coast  at 
1.5  km  above  the  Adriatic  Sea. 

Coupled  and  uncoupled  simulations  of  the  bora  were 
performed  using  COAMPS  with  a  horizontal  grid  in¬ 
crement  of  5  km.  The  strong  downslope  winds  in  the 
simulation  extended  westward  over  the  Adriatic  and  re¬ 
sulted  in  a  significant  enhancement  to  the  surface  rough¬ 
ness  as  a  result  of  the  growing  waves  in  the  fetch-limited 
conditions,  as  shown  in  Fig.  12a.  Plumes  of  large  surface 
roughness  extend  across  the  Adriatic  and  are  associated 
with  the  fetch-limited  downslope  flow.  However,  in  spite 
of  the  larger  roughness  used  in  the  atmospheric  model 
due  to  the  coupled  feedback,  in  general,  the  simulated 
atmospheric  and  wave  fields  are  quite  comparable.  For 


example,  the  wind  speed  difference  at  10  m  between 
the  coupled  and  uncoupled  simulations,  shown  in  Fig. 
12b,  differs  by  less  than  3.0  m  s-1,  which  is  compar¬ 
atively  small  in  relation  to  the  maximum  bora  speed  of 
25  m  s  !  near  the  surface.  Overall  the  wave  coupling 
appears  to  have  minimal  impact  on  the  atmospheric 
wind  and  thermal  fields.  The  timescale  in  this  case  is 
likely  too  short  to  produce  a  significant  two-way  inter¬ 
action  response  between  the  atmosphere  and  ocean  sur¬ 
face  waves.  This  result  is  in  overall  qualitative  agree¬ 
ment  with  Janssen  and  Viterbo  (1996)  who  found  that 
the  impact  of  wave-atmosphere  coupling  is  significant 
for  longer  timescale  processes,  such  as  those  that  affect 
the  atmospheric  and  wave  climate. 

4.  Summary  and  conclusions 

In  this  study,  three  tropical  cyclones  and  a  topograph¬ 
ically  forced  bora  event  over  the  Mediterranean  Sea 
were  simulated  using  a  coupled  atmosphere-ocean  wave 
model,  COAMPS.  The  model  consists  of  a  nonhydro¬ 
static  atmospheric  model,  COAMPS,  coupled  with  a 
wave  model,  WAM.  The  coupled  modeling  system  is 
capable  of  enhancing  surface  stress,  due  to  young  ocean 


1  .0 

0.5 

0.0 


Fig.  12.  Simulated  (a)  roughness  length  difference  (cm)  and  (b)  10-m  wind  speed  difference  (ms  ')  between  the  coupled  and  uncoupled 

simulations  valid  at  1200  UTC  7  Nov  1999  (12  h). 
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waves  in  fetch-limited  conditions,  beyond  values  yield¬ 
ed  by  the  conventional  Charnock  scaling  relationship. 
The  models  were  integrated  simultaneously  and  on  iden¬ 
tical  grids  with  horizontal  grid  increments  ranging  from 
12  to  5  km.  The  impact  of  the  wave-induced  stress  was 
largest  for  the  strongest  of  the  three  tropical  cyclones, 
Mitch.  The  central  pressure  in  the  coupled  simulation 
in  this  particular  case  was  8  hPa  more  intense  than  the 
uncoupled  simulation  due  to  larger  surface  heat  and 
moisture  fluxes  that  significantly  augmented  the  precip¬ 
itation.  This  intensification  apparently  resulted  primarily 
from  a  weak  dependence  of  heat  and  moisture  flux  on 
the  wave  state.  Other  tropical  cyclone  simulations  in¬ 
dicated  less  sensitivity  to  wave-induced  effects.  In  gen¬ 
eral,  the  low-level  wind  speeds  were  insensitive  to  the 
coupling  in  spite  of  a  Charnock  parameter  increase  by 
more  than  an  order  of  magnitude.  Maximum  low-level 
wind  speeds  typically  were  2-3  m  s-1  less  due  to  the 
feedback  of  ocean  wave-induced  stress.  However,  local 
differences  in  excess  of  10ms_1  were  found.  Significant 
wave  height  maxima  were  reduced  by  approximately 
10%  in  the  coupled  simulations  due  to  the  enhanced 
roughness  associated  with  the  young  ocean  waves.  The 
reduced  significant  wave  heights  simulated  by  the  cou¬ 
pled  model  were  in  closer  agreement  with  the  SRA  ob¬ 
servations  for  Tropical  Cyclone  Bonnie  than  for  the  un¬ 
coupled  simulation.  The  simulated  tropical  cyclone 
tracks  were  insensitive  to  wave  coupling  effects.  The 
coupled  model  indicated  reasonable  skill  in  replicating 
some  complex  wave  spectra  characteristics,  such  as  mul¬ 
timodal  structures,  as  observed  by  SRA  for  Tropical 
Cyclone  Bonnie.  Topographically  forced  flows  in  the 
coastal  zone,  such  as  the  bora,  are  characterized  by  oc¬ 
casionally  high  near-surface  wind  speeds  and  fetch-lim¬ 
ited  conditions.  However,  in  this  study,  the  bora  event 
was  the  least  sensitive  to  the  ocean-wave-enhanced 
roughness  of  all  the  cases  considered.  For  example,  in 
spite  of  a  substantial  increase  in  the  surface  roughness 
associated  with  the  young  ocean  waves,  the  near-surface 
wind  speeds  differ  by  less  than  3  ms-1  between  the 
uncoupled  and  coupled  simulations. 

It  has  become  increasingly  apparent  that  surface 
ocean  wave  effects  due  to  wave  age  factors  is  a  key 
component  of  the  air-sea  exchange.  Incorporation  of 
surface  ocean  wave  effects  into  numerical  models  may 
be  necessary  to  achieve  accurate  surface  fluxes,  partic¬ 
ularly  in  high  wind  conditions.  A  number  of  issues  re¬ 
main  unresolved  with  respect  to  air-sea  interaction  in 
high  wind  conditions.  The  accuracy  of  wind  energy 
transfer,  dissipation  and  nonlinearity  representations  in 
wave  models,  surface  flux  exchange  at  the  air-sea  in¬ 
terface  in  atmospheric  models,  and  air-ocean  coupling 
methodologies  need  to  be  assessed.  In  the  high  wind 
regime  in  particular,  gustiness  effects  and  mesoscale 
transients  need  to  be  considered  in  order  to  improve 
surface  wind  stress  forcing.  Finally,  an  improved  rep¬ 
resentation  of  the  bulk  effect  of  sea  spray  on  the  fluxes 
of  heat  and  moisture  may  be  required  before  realistic 


simulations  of  air-sea  interaction  at  high  wind  speeds 
can  be  achieved. 
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